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CHAPTER  I 
INTRODUCTION 

Elasticity 

In  the  study  of  elastic  properties  of  materials  forming  regions 
vdth  greater  than  microscopic  dimensions,  the  molecular  nature  of  the 
material  Is  disregarded  and  replaced  by  a  continuous  mathematical 
model.  The  problem  does  not  deal  with  movements  In  the  body  as  a 
whole,  but  in  relative  movement  of  points  within  the  body,  vrtilch 
movements  are  called  strain. 

The  forces  that  produce  strain  are  referred  to  as  stresses  and 
are  divided  Into  two  classes.  Normal  stresses  tending  to  elongate  or 
shorten  a  segment  are  called  tensile  or  compressive  stresses,  respec- 
tively. A  stress  tending  to  deform  an  angle  within  the  body  Is 
called  a  shear  stress. 

To  describe  the  stresses  associated  with  rectsingular  coordinates, 
three  sjmibols,  Gx»  CTy»  and  (j^,   are  necessary  for  the  normal  streses; 
and  three  symbols,  f^u,      1x2*   and  iyg,   for  the  shearing  stresses. 
Similarly,  for  the  strains,  we  use  Cx,  Cy,  C^,  ffxy  >      lf%z>  ^^^     VY^' 

Moreover,  it  may  be  shown  that  the  stress  coordinates  are  co- 
variant  tensors  of  the  second  kind,  as  are  also  the  strain  coordi- 
nates (6,  p.  275).  This  fact  simplifies  the  changing  of  spatial 
coordinates. 


To  the  theory  Hooke  (8,  p.  l)  contributed  his  hypothesis  that 
there  exists  a  linear  relationship  between  stresses  and  strains  when 
they  are  of  infinitesimal  order.  Below  are  the  relations  of  Hooke *s 
Law  as  they  are  regularly  now  given. 

(1.1)        €y  =  -^[cTy  -Via,  +  crjj ,  j,^=± r„  , 

where  E,  G,  and  V   are  elastic  constants  pertaining  to  the  material 
in  question,  and  are  called  the  modulus  of  elasticity  in  tension,  the 
modulus  of  elasticity  in  shear,  and  Poisson's  ratio,  respectively. 

Equilibrium  within  the  elastic  body  demands  that  the  following 
equations,  called  the  differential  equations  of  equilibrium  (9,  p.  229), 
be  satisfied: 

9x   dy  dz  ^  ^ 

where  X,  Y,  and  Z  are  components  of  any  body  forces  present. 

Moreover,  relations  between  displacements  and  strains  reqiiire  that 

<3y^     5x^    <^x9y 


<  3a 


(1.3) 


SySi.  3V  I     o)x  c^y  c^z 

<3xdz         <9y  I    ^x         5y  9^ 

^  >9x9y        ^il   ax         Qy         9i 

which  equations  are  called  the  conditions  of  compatibility  (8,  p.  28). 

To  solve  a  general  problem  in  three  dimensions  when  a  body  is 
submitted  to  the  action  of  given  forces,  it  is  necessary  to  solve 
equation  (1.2),  subject  to  boundary  conditions,  and  to  obtain  a  set 
of  stresses  which  produce  a  corresponding  set  of  strains  that  satisfy 
equation  (1.3).  Anyone  familiar  with  such  problems  will  realize  that 
this  is  a  tremenduously  complicated  and  difficult  task,  except  in  the 
case  of  the  simplest  stress  distributions.  One  must  also  realize  that 
it  is  here  assumed  that  the  material  has  the  same  elastic  properties 
in  each  direction.  This  means  that  a  completely  general  solution  in 
elasticity  is  difficult  to  the  point  of  impossibility.  It  is,  there- 
fore, customary  to  introduce  certain  simplifying  assunqstions  that  re- 
rove  many  conditions  on  the  general  solution  and  still  yield  results 
that  approximate  many  situations  closely  enough  to  be  very  useful. 
In  the  remainder  of  this  work,  the  subject  of  plane-elasticity  will 
be  followed. 

Let  it  he  assumed  that  the  dimension  of  the  body  in  the  z  direc- 
tion is  small  relative  to  the  other  two.  Suppose,  also,  that  the 


loading  on  the  edges  of  such  a  body  is  distributed  uniformly  along  the 
thickness  and  that  its  faces  are  not  loaded.  A  large  plate  of  steel, 
loaded  on  the  edges,  wovild  be  an  example  of  this  problem.  With  this 
loading  pattern,  O^,  Txj,  and  Tyj  are  zero  on  the  boundary,  and  it  is 
also  assumed  that  they  are  zero  inside  the  body.  Likewise,  we  assvune 
that  the  remaining  stresses  are  functions  of  x  and  y,  but  not  of  z. 
Such  a  situation  is  called  plane  stress. 

Under  these  assumptions,  it  is  seen  that  the  equations  (l.l), 
(1.2),  and  (1.3)  reduce  to 


(1.1a) 


(1.2a) 


dX     d  y 


Another  similar  simple  case  called  plain  strain  is  derived  by 
assuming  that  the  z-direction  is  very  large  with  respect  to  the  other 
two.  Consider  the  case  when  forces  load  the  body  perpendicular  to  the 
z-direction  and  are  uniformly  distributed  along  its  length.  Then,  it 
is  assumed  that  there  is  no  displacement  in  the  z-direction;  hence, 
there  are  no  strains  in  this  direction.  As  before,  the  stresses  and 
strains  will  depend  upon  x  and  y  only.  Since  C2,  ^^,  and  f^^  ^^® 


zero,  then  (l.l)  reduces  to 

e.=  i[cr.  -vicryl <y^]  ^ 

(l.lb) 

showing  that  C^  may  be  expressed  in  terms  of  (5^  and  d^.   Eqxiatione 
(l.lb)  then  become 

Equations  (1.2)  and  (1.3)  reduce  to  precisely  the  same  form  as  (l.2a) 
and  (l.3a)  above. 

It  can  here  be  seen  that  plane  stress  and  plane  strain  are  iden- 
tical mathematically,  and  differ  only  in  the  elastic  constants  enqjloyed 
in  equations  (l.la)  and  (l.lc).  These  two  can  then  be  grouped  to- 
gether under  the  subject  of  two-dimensional  problems. 


Two-Dlmenslonal  Problems 

As  has  been  seen,  in  order  to  solve  a  two-dimensional  problem, 
we  must  satisfy  at  least  three  sets  of  conditions.  These  are  the 
equilibrium  equations  (l.2a),  the  conditions  of  compatibility  (l.3a), 
and  the  boundary  conditions. 

G.  B.  Airy  (9,  p.  26)  first  introduced  the  stress  function  (p 
such  that,  when  there  are  no  body  forces, 

<-'       1f^-^'^  S^-^r^         S$=-^«/. 

This  function  satisfies  the  eqiiilibrium  conditions. 

J.  C.  Maxwell  (7,  p.  105)  found  that  if  the  function  0  satis- 
fies the  bi-harmonic  equation,  then  the  condition  of  compatibility 
is  also  satisfied.  The  problem  now  resolves  into  finding  a  (p   satis- 
fying 

(1.5)  1:^  +  2^+^  =  0, 

and  also  satisfying  the  boundary  conditions. 

It  is  found  that  a  general  solution  of  the  bi-harmonic  eq\iations 
is 

This  may  be  written 

(1.0  0=  R[zV^fe)+x(z)]  ^ 

where  z  =  x  +  iy  and  y   and  X   are  complex  functions  of  z,  analytic  in 
the  region  under  consideration.  Some  elementary  manipulations  lead  us 
to  two  relations  that  are  used  to  evaluate  the  stresses,  avoiding  the 


necessity  of  finding  the  real  part  of  Cf),     They  are 

(1.7)       CJ,  +  a   =  2(y/'*F)  , 

(1.8)  cr^- o«+2;t;,  =  2(2V^"+x"). 

It  is  also  found  that  the  resultant  force  acting  along  a  boundary, 
vfhere  X  and  Y  are  the  restiltants  in  the  x  and  y  directions,  respec- 
tively, and  A  and  B  are  the  ends  of  the  boiindary,  is  (7,  p.  IIA) 


(1.9)       x-iY  =  -iLv^-^zy^'+x']" 


It  follows  that,  along  an  unloaded  bovindary,  this  quantity  must 
be  constant.  Written  in  the  way  that  we  shall  repeatedly  use,  it  is 

(1.10)  [y^  +  H'^-fX'']    =     H  , 

where  H  is  a  complex  constant. 

With  the  foregoing  relations,  (1.7)  through  (l.lO),  certain  inter- 
secting boundary  problems  will  be  attacked.  It  is  most  important,  how- 
ever,  to  discuss  the  functions  ^f  and  X ,  the  complex  potentials, 
before  proceeding. 

The  Complex  Potentials 
Several  rather  general  statements  may  be  made  about  llf  and  X. 
It  has  been  noted  that  they  must  be  analytic  in  the  region  under  con- 
sideration. They  may  have  singtaarities  in  deleted  regions,  and  there- 
fore, may  be  expanded  in  a  Laurent  series  about  these  regions,  except 
for  any  terms  that  might  produce  multi-valued  stresses  or  displace- 
ments. It  is  to  be  noted  that  the  imaginary  part  of  the  function 


[glZ/'  +  xl  may  be  mult i -valued.  Muskhelishvili  (7,  p.  121)  shows  that 
XJf  and  X  iMst  have  the  following  forms : 

(1.11)  ^  =    £a,Loc3U-2J    ^     f\ 


m 


(1.12)  X    =  5]6,Lo<^U-z,)  *   2X!^KLoQ(i-z,)  +  X*, 

K»l  K=:l 

where  Aj^,  B,  ,  and  C.  are  con^jlex  constants,  and  Uf   and  X  ai"e  holo- 
morphic  functions.  Zj^  is  a  point  in  the  k^^  deleted  neighborhood. 

Since  the  stresses  do  not  depend  on  X  but  on  X  and  higher  de- 
rivatives, the  X  form  is  here  inserted. 

m 
(1.13)        X'  =  ^CkLoc.(2-zJ  4-  X,  , 

where  Xi  is  holomorphlc  in  the  region  of  the  problem. 

I     It  It 

Since  Uf  and  X,  are  holomorphic,  they  may  be  freely  expanded  in 
laurent  expansions  about  points  in  deleted  neighborhoods. 


CHAPTER  II 
THE  MOTIVATION  FOR  THE  METHOD 

As  can  be  readily  seen,  the  complex  functions  \U   and  X.   yield 
results  that  are  sometimeB  readily  handled  by  the  method  of  undeter- 
mined coefficients.  The  usual  procedure  is  to  insert  the  boundary 
conditions  in  an  appropriate  relation,  such  as  (l.lO),  and  to  use 
the  fact  that  t>ro  expressions  are  identical  if  coefficients  of  like 
powers  of  z  are  equal.  This  gives  a  series  of  equations  which  are 
to  be  solved  for  the  vinknown  coefficients. 

In  the  past,  attempts  to  solve  certain  problems  by  the  method  of 
undetermined  coefficients  have  failed  to  yield  resiilts.  One  of  these 
cases,  the  semi -circular  notch  problem,  which  has  been  solved  using 
the  methods  of  real  Arariables  (5),  and  later  generalized  by  the  same 
methods  (3)»  will  now  be  discussed  to  show  where  the  difficulty  arises. 
As  shall  be  seen,  inconsistent  coefficient  equations  result  from  the 
presence  of  a  multiplicity  of  boundaries. 

The  notch  in  question  has  been  ^^^^ 

oriented,  as  at  the  right,  so  that 

P 

it  becomes,  by  appropriate  choice  of   -< 


y.j 


p 

— ^ 


axes  and  scale,  the  half  of  the  unit  „ 

>if    7 

circle  for  which  y  is  positive.  This  Figure  1 

portion  of  the  unit  circle  will  be  called  the  half  unit  circle. 
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The  X-axis  is  the  edge  of  our  semi -infinite  sheet.  Since  this 
can  be  thought  of  as  an  infinite  sheet  with  a  circular  hole  coupled 
vdth  some  distributed  force  superimposed  on  the  x-axis,  we  can  ex- 
pect the  complex  potentials  to  have  singularities  only  within  the 
unit  circle.  Then  Jt/ and  X    will  asstune  the  following  forms: 


Moreover,  since  the  stresses  are  bounded  at  infinity,  we  can  immedi- 
ately eliminate  all  positive  powers  of  z  except  the  first.  Then,  also 
deleting  the  constant  terms  not  affecting  stresses,  we  have 


oo 


Xj/     =  ALo<52   +  az   +  ^Qn2*''  , 

oo 

^  r=l 

or,  in  a  better  form  for  our  purposes, 

oo 

(2.1)  ^   -  ALo^z  1-  (A'-^A"0^   -^'Y^a^z-''  , 

OO 

(2.2)  X'=    BL03E    +(E)'^B")£     +Eki"    . 

The  uniform  tension  P  in  the  x-direction  fixes  the  coefficients 
of  z,  for  from  (1.5) 

(a.^cs,)       =  ^{\ff'-^y7')    =  4A'  =  P   , 

since  Ol  =  P   and  c^y  =  O  .  Then  A'  =  ^« 

*  \oa  /loo  '♦• 


n 


Also,   from  (1.6) 


oo 


80  that     B'=-y,     B"  =  O  ,   since    Xcy 


=  zi^xl;"^  X"]--  2(B'-iB")  =  -P  , 

=  0.  A"  is  set  equal  to  zero. 


since  it  does  not  affect  the  stresses  and  since  the  movement  about  the 
origin  also  vanishes. 

The  final  form  of  \u   and  X   to  be  used  in  succeeding  work  is 


(2.1a) 

(2.2a)  X'   -   BLoc^^  -!:£  +^^,2"^. 


oo 


n=i 


Since  no  force  acts  on  the  x-axis,  then  (l.lO)  holds  and 


=    H. 


(1.10a)  [!//+Z)//'^X'] 

Substituting  from  (2.1a)  and  (2.2a)  yields 

oo  T-  **  \       _ 


(2.3) 


-^-,^^k,r"  =  H 


tn  =  l 


The  above  expression  must  be  equal  to  a  constant  along  the  line  z  =  x. 
However,  on  this  line,  z  =  z  =  x,  and  the  foregoing  becomes 

(A  +  B)Lo<^x  -(-a  +^(Q„-nQ^+b^)x"''  =H. 

This  is  satisfied  if   B  -  "A  and   D^=  na^-a^^  or 

B  =  - A  , 

by  taking  conjugates  of  both  sides.  Then  (2.1a)  and  (2.2a)  simplify  to 


12 
(2.1b)  )//'  =ALo3£^  ^z  +^Q.2-"   , 

h  =  l 
go 

(2.2b)  X'    =-ALp<^i-|£  +2](nQ„-a^)i~". 


n=l 


Now,  considering  the  unit  circle,  it  is  seen  that  z  =  e    along 
this  curve,  since  the  modulus  R  is  equal  to  one.  For  simplicity, 
e^®  will  be  called  Cf*      a  =  e'^  =  Cf"'  ,  so  that  3  -    O"'  .  In- 
serting this  boundary  in  (2.3)  leads  to 

H=  ALoocr-AL^cr-' +|o'-|a'-'  -^Acr^ 

The  usual  attack  would  yield  the  following  equations  when  coeffi- 
cients are  set  equal: 

(a)  A  +  A  =  0  (d)  A  +  b2  =  0 

(b)  f  *  b^  =  0  (e)   a^  =  0,     n  >  1 

(c)  a^  -  2  =  0  (f)     bj^  -  (n-z)a^_2  =  0  ,         n>2 
Careful  scnitir^  of  these  eqxiations  shows  them  to  be  inconsistent 

when  considered  with  equations  (2./*).  Note  that  b^  =  na^  -  a^  coupled 
with  (f)  above  shows  a  contradiction,  in  that,  from  the  former  b^  =  0 
since  it  depends  on  ao,  whereas,  the  latter  yields  b„  as  dependent  on 
a^  and  therefore  bo  ?^  0.  There  is  agreement  at  this  point,  however, 
that  A  must  vanish. 


13 

The  usual  method  has  therefore  failed  to  produce  fruitful  results. 
This,  however,  poses  a  more  far-reaching  problem:  that  is,  the  exist- 
ence of  the  problem,  at  least  theoretically,  would  show  that  the  very 
fujictions,  which  have  led  into  this  inconsistency,  must  in  some  way 
give  a  solution.  This  would  indicate  that,  somewhere  in  the  process, 
there  must  have  been  some  step  taken  that  is  based  on  the  grounds  of 
sufficiency  but  not  necessity.  That  step  is  discussed  in  the  next 
chapter. 

In  summary,  there  is  need  to  investigate  the  methods  commonly 
used  because  they  lead  to  no  solution,  where  it  can  demonstrated  that 
a  solution  exists.  This  investigation  will  begin  with  a  look  at 
identities  on  the  half  unit  circle. 


CHAPTER  III 
IDENTITIES  ON  THE  HALF  UNIT  CIRCLE 

The  half  unit  circle  vd.ll,  as  has  already  been  noted,  in  the  fol- 
lovring  context,  refer  to  the  upper  half  of  a  unit  circle  centered  at 

the  origin,  that  is  the  part  of  the  unit  circle  with  positive  ordinate. 

i0 
If  z  is  expressed  as  Re  ,  then  on  this  portion  of  the  unit  circle 

R  =  1   , 

0  <  e<  tt  . 

But  Euler»s  faniiliar  formula  states  that 

16 
e   =  cos  9  +  isin  9, 

or,  by  rearrangements. 


;e  _     10 
(3.1)        sine  =  -^   ^.  ^ 


COS&  =  -^ — ^-^ —  . 

'ft 

Previously,  the  symbol  <?=  e'   was  introduced,  and  it  was  there 
noted  that  0-'®=5=cr''.  Also  e"'®  =  C^ ,     Then 


Sin(mQ)  = 


o""-  d" 


(3.2) 


•> 


C05(mQ)  = 


14 


15 

In  the  region  under  consideration  sin  dd   and  cos  m©,  where  m  is 
a  positive  integer,  can  be  expanded  by  an  ordinary  Fourier  series  into 
an  expression  involving  only  cos  nG  and  sin  nB,  respectively.  This 
is  not  true  of  the  whole  circle  where  uniqueness  demands  that  the  ex- 
pansions for  sin  ni8  and  cos  mG  be  precisely  sin  mS  and  cos  mQ,  respec- 
tively. 

By  appropriate  Fourier  theorems  (l,  pp.  57,  70),  we  obtain  ex- 
pansions of  the  form 

cr"'-o^"    -  -^ -h    fc-^-t-^^   ,      o<e<TT, 


n=i 


where 


Co  =     |/    Sin(m9)d0      , 


'O 

'IT 

2 


and 


Cn=    i        Sin(me)co5(na)d0  ,     n=  1,2,3,-, 


n^l 


where 


D^   =  |:/cos^e)5in(ne)cie5     n  =  i,2,3,---. 

•'a 


It  is  found  by  simple  integration  that,  when  m  =  2k  -  1,  then 

^     n=  0,1,2,3,..., 


r  -  gf  2(gK-i) 


C^n-P  O 


16 


and,  when  m  =  2k,  then 


^2n-.   1T  L(^K)^-Un-I)2j  ,   n=i,<;,J,' 


Also  if  m  =  2k  -  1,  then 


and,  when  m  =  2k,  then 


D^n  =  0  , 


f.    _  2  r   £C2n-l) 

^^^-^  -rr  L(2n-i)2-UK> 


■ji   n=  1,2,3, - 


These  lead  to  the  following  relations: 


i.K~\         =7  tK- 


cr^"^-'-^ 


2i 

c 

a; 

a 

e,K-i     •=>2»<>-t 

a 

0 

2K           r=  2K 

1    0-2n  ,  »2.>i 


n=r 


a'^^-Ha 


=  J-  V  r   a  UK)   1 
_  a  yf    2(gn)    1 


2(2^.-1) 


L(dn-I)^-(^K)^ 


^a.-i 

+  5^^'^-' 

2 

a^" 

-a" 

2 

^2... 

_^2n-. 

2i 


where  0  <  6  <  IT  . 
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By  elementary  manipulation,  we  may  obtain 


(a) 


o        -  O  -    ^ 


(c5*"+^^» 


(3.3) 


nsl 


(c)     cr     H-C7       =  irLaK-i)2-(^n^^^^     -cr   j, 

oo 

(d)        (^        -hc5  =lf  L  (ZK)--(2Ln-l)i^'^  ^        ^- 


hTi 


Properties  of  infinite  series  would  allow  series  to  be  added, 
with  the  assurance  that  the  resulting  series  will  also  converge  in  the 
same  region.  If  we  add  (a)  and  (c)  from  above,  we  obtain 


cr^"-'=-|J-< 


=^^ 


ZK- 


+ 


2K-I 


■I 

r\=i 

■E 


2r. 


-h 


1' 


By  Abel»B  test   (4,   p.   SO), 


Consider  the  series  ^  -  raK-i^  -[Zr\\     ^^' 
this  series  will  converge  at  all  points  on  the  unit  circle  except 
possibly  at  z  =  1,  since  the  coefficients  of  the  negative  of  this 
series  form,  after  the  2k^^  term,  a  decreasing  sequence  of  positive 
terms  approaching  zero  as  n  increases  without  limit.  A  similar 
statement  may  be  made  about  the  second  series  given  in  (3.2). 

Since  both  series  converge  on  the  unit  circle  where  O  <  6  <  rr  , 
the  terms  may  be  separated  as  long  as  use  of  the  series  is  confined  to 
this  region.  Then 


^"-'-f 


S-2M 


2K-I 


+y tf'"    +y  — ^ 

^  Z_.  2K-l-^»n    Z-.   2K-I4- 


Zty 


0  < 


e  <'TT 
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If  we  use  c:5=  C3'"'   and  substitute  -n  for  n  In  the  second  stun- 
mation,  we  have 


^aK-,  _  2; 


■  Zn 


-  00 


tf  I  2K- 


ZM 


Since 


iK-  I 


(2K-I)  -  (en) 


^zL  2K-I  -2n  ^Z^  2K-l-^n 
n=i  n=-l 

,  summations  may  be  combined  giving 


B  =  0 


■  an 


(3.4a) 

Now,  adding  (b)  and  (d)  of  (3.1)  produces 


n=l 


'  iT  Z-L  2K-l2n-l)  ^  2K+(2.n-l)J 

As  before,  the  two  series,  made  up  of  the  first  and  second  terms 
of  each  term  of  this  series,  each  converge  in  the  region  0  <  S-  <  Tf 
Then  the  same  rearrangement  yields 


,2K 


=fE- 


ZK-Un-I)  "^ 


y : 


:='2n-l 


lETT^iy 


Substituting  in  the  second  summation  -(2n  -  l)  for  2n  -  1  gives 


^2K 


n=i  r»=o 


and  combining  summations  obtains 


&n-\ 


(3.5a)       Cf  =-|^E  2K-V-0  • 

From  these  equations,  (3.4a)  and  (3.5a),  it  is  seen  that  there 
exists  a  relation  between  power  terms  in  Cf   on  the  unit  circle. 


19 

Since  this  is  the  case,  then  there  might  be  several  vraiys  for  a  function 
to  vanish  on  the  half  unit  circle.  It  might  vanish  because  the  coef- 
ficient of  each  power  vanishes,  by  certain  infinite  sums  of  coefficients 
vanishing  through  the  relation  (3.4a)  and  (3.5a),  or  through  some 
other  relation  between  one  power  of  CS   and  others  in  an  isolated  re- 
gion. The  second  possibility  will  be  used  in  the  sequel. 


CMPTER  IV 
THE  METHOD 

On  the  half  unit  circle,  it  has  been  shown  that  there  exists  a 
relationship  between  powers  of  d".  These  relations  are 

a-  _  2;  V^  CT^" 


(2K-0-2r> 


n  =  - 


(3.5a)  '^     "  ir /L  UK)-Uh-i)      o<e<Tr. 


nx-00 


The  form  of  the  above  equations  would  imply  that  the  imaginary- 
part  of  the  coefficient  of  odd  powers  affects  the  real  part  of  the  coef- 
ficients of  even  powers,  and,  similarly,  that  the  imaginary  part  of 
the  coefficient  of  even  powers  affects  the  real  part  of  the  coefficients 
of  odd  powers.  This  would  suggest  dividing  the  expressions  which  are 
to  be  evaluated  into  three  parts,  i.e.,  powers  of  C^  with  real  coef- 
ficients, even  powers  of  (3   with  imaginary  coefficients,  and  odd  powers 
of  CX  with  imaginary  coefficients.  Then  (3.4a)  and  (3.5a)  can  be  ap- 
plied, the  imaginary  unit  therein  causing  all  coefficients  involved  to 
become  real.  Gathering  the  coefficients  of  like  powers  will  then  pro- 
duce real  coefficient  expressions  which  are  usually  infinite  series. 
Then,  setting  each  such  expression  equal  to  zero,  or  to  corresponding 
boundary  conditions,  will  remove  all  stresses,  or  otherwise  match 
given  boundary  conditions. 
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To  illustrate  the  use  of  (3.4a)  and  (3.5a),  application  vd.ll  be 
made  to  the  circular  notch  problem  which  was  started  in  the  second 
chapter  of  this  vrork.  Since  this  problem  has  been  solved,  this  ap- 
plication vfill  serve  to  corroborate  the  method,  as  well  as  to  increase 
the  evidence  of  the  accuracy  of  the  first  solution. 

This  method  is  applicable  to  any  two-boundary  problem  in  which 
one  boundary  is  mapped  onto  the  half  unit  circle.  If  the  boundaries 
are  disjoint,  then  application  becomes  very  difficult.  The  easiest 
and  most  straightforward  application  is  to  boundaries  which  intersect, 
such  that  one  can  be  mapped  on  the  semi-circle  and  the  other  on  the 
X-axis.  Such  is  the  problem  of  Chapter  VII,  which  heretofore  has 
been  unsolved. 

It  should  be  here  added  that,  although  the  results  are  ueiially 
infinite  series  in  undetermined  coefficients,  they  tend  to  assume  sim- 
ple forms  which  possibly  could  be  treated  generally  to  get  exact  re- 
stilts,  or,  at  least,  to  corroborate  hypotheses  concerning  general 
solutions. 


CHAPTER  V 
THE  SEMI -CIRCULAR  NOTCH 

The  method  just  described  will  now  be  applied  to  the  circular 

notch  problem  which  was  begun  in  Chapter  II.  The  complex  potentials 

were  restricted  in  form,  by  the  vanishing  of  the  logarithmic  term 

and  the  tension  at  infinity,  to 

oo 


f  -^i+Zo„^-\ 


(5.1)  oo 

To  remove  stress  on  the  x-axis,  it  was  necessary  that 

(2.4)       b  =  non  -  Qn  . 

The  boundary  conditions  must  now  be  satisfied  on  the  semi -circular 
arc.  First,  substituting  (2.4)  in  (5.1)  yields 


(5.1a) 


Then 


h=l 


(1.10)  y/  -f  zifr'  +  x' 

where  H  is  a  complex  constant. 


=  H  , 

2=^ 
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Inserting  (5.1a)  in  (l.lO)  gives 


Z  =  cr 


nil  ^         •^sl 


n=l 


Now,  if  we  let  a  =  c  +  id  ,  then 

n    n     n 


hsi 


=  -r^  •    z. 


p  «-i 


oo 


nn 


n=i 


This  separation  is  permissible  since  the  convergence  of  a  complex 
series  insures  the  convergence  of  the  series  composed  of  its  real  and 
imaginary  coefficients,  respectively. 

In  the  eqimtion  at  the  top  of  the  next  page,  the  coefficients  of 
some  of  the  powers  of  tf  are  gathered,  and  some  of  the  summations  are 
altered  and  combined. 
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Letting  the  part  of  the  above  equation  with  real  coefficient  be  called 
F{ci),   and  dividing  the  imaginary  coefficient  part  into  even  and  odd 
powers  of  C?  gives 


n=l 


r\z.Z 
oo 


The  preceeding  expression  is  now  in  convenient  form  to  apply 


formulas  (3 •4a)  and  (3 •5a)  to  obtain 


.2K-I 


ZK-I 


oo 


K=-oo 


■D^^-(^)E-^;;?5^   -^.d,!^)! 


ZK 


oo 


K=-<X5 


K=-oo 


I  -  2K 


EK-I 


n  =  *  K--eo 

«o        oo        2K  ""^       EK 


n=i 


K.C-00 
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The  expression  now  has  all  real  coefficients.  By  referring  back 
to  the  value  of  F{0), 


OO  oo 

^2  r""  ^-*i  ^TT 


it  is  possible  to  pick  out  coefficients  of  various  powers  of  C  ,  and 
to  set  each  one  equal  to  zero,  since  the  left  member  of  the  eqiiation, 
H,  is  a  complex  constant. 

The  coefficient  of  the  first  power  of  CT   is 


_P_ 

2. 


-»- 


MZ-^^  -n(^-')d.„-(^"-^)d..-J^ 


•rial 


P:E 


3di 


2- 


_     _P_    .    Z_\y     dan        ,Y"Un-el)d.n     _y  (gn)d2n 


JtsL 


z~  I 


-5-  + 


in 


■1 

oo 


1  Er.-! 


2n        ]J 


a 


.      2    J    8    J        ,V       ^(^^"^       A      ' 


(5.2)         _     -^    +7rLten^^-|      ^^" 


=    0   . 


26 


The  coefficients  of  d^  and  all  higher  odd  powers  take  the  form 
for  the   (2k  -  l)^^  power: 


0  =(2K-2)C,,.,-(2K-3)q,.3  +  ^. 


n:2 


3d 


3dz         .V^ d, 

-Uk-1)      Z^  2n  +  ( 


IE 


+yr ^ ^n-l         1  d  V       ^^  dm 


n»i 


''^^-^'^--■-^^^-^'^---'^^Ic.-um'jI^-uk-oj 


+ 


OO 

n«2 


4Un)(£r,  +  0d^ 


C(2.»^)^-UK-i)^JU«+a-C2K-0] 


(5.3) 


OK} 


The  coefficient  of   CJ"*  is 

'^  '2n-l       Z_i  2n+l  2+  I 


a    ■  ^'■''  TT   1^2, 


—    £.4.r  4.  2 


_      P   ,  r   ,   a  f  8d^     ■  V        4(ah)dan 


(5.A)  =   -£■  +  '-'+  Tr|L   (En-l)(an-3)/  • 
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The  coefficient  of  O" 


2K-I 


,  when  k  >  1,  is 


0   =  C^K-,-*-  tT 


3d, 


Ln=l 


C2XHI 


n=2 


(2n-H)den-(2r^-2)da„-, 

2h  +2K-I 


S+Z.K-1 


=  c 


2K-I 


IT  jLz  +  UK-l)  "*■  2-C2.K-n  ~  4  +  (2K-l)J  ' 


nl 2njtl__1J    V    £ndtn 


n-Z. 


=  c 


iK-1 


+  ^1 


2+dt 


n»2 

oo 


(2K-I) 
4(^h)(2h-H)d 


ic_ 


^-(^KW)^J[4f2K-(|  '*'^4l[W-(2>^-l)]fe^^^+^'^-fl 


r*;  «;^    —  r         J-  2  V  4(£n)(£n-H)d2n       , 

K=2,3,4,.--. 

It  is  seen  that  the  eqiiations  (5.2)  through  (5.5)  contain  only 
even  numbered  d»8  and  odd  numbered  c»s.  The  coefficient  equations 
obtained  from  the  even  powers  of  d  contain  only  odd  numbered  d»s  and 
even  numbered  c*8.  Furthermore,  these  equations  contain  no  constants, 
and  therefore,  form  a  homogeneous  system  of  equations.  On  this 
basis,  we  may  solve  all  that  system  merely  be  lotting  all  odd  d  and 
even  c  coefficients  be  equal  to  zero.  The  solutions  thus  obtained 
agrees  with  the  solution  of  this  problem  by  Maunsell  (5)  as  is 
shown  in  Chapter  VII. 

A  system  of  forty  equations  in  forty  of  these  coefficients  was 
composed  in  the  following  order: 

1^''^  and  2^  equations  from  (5*4)  and  (5.2)  respectively, 


f^,   5th  ^  y'^^ ,  •  •  • ,  39*-^  from  equation  (5.5), 


.th 


k     ,   6^^,   8^^,"*,  40^^  from  equation  (5.3)  . 
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This  order  produces  a  system  with  the  coefficient  matrix  diagonal 
that  is  large  compared  to  the  rest  of  the  matrix.  Its  simplicity  of 
form  allows  rapid  machine  computation,  and  the  heavy  diagonal  greatly 
improves  accuracy.  This  system  was  solved  to  eight  significant  digits 
on  an  IBM  650  calc\ilator  in  the  Statistical  Laboratory  of  the  Univer- 
sity of  Florida.  Below  are  the  first  thirty  rounded  coefficients  of 
y/"  and  X'. 

TABLE  I 
COEFFICIENTS  0?  ^jf  MD  X  ' 

n  yr  X' 

1 +0.7924848  +0.0000000 

2 -  .33141211  -  .99423631 

3 -  .2173459  -  .4346918 

4 +  .08476891  +  .42384461 

5 +  .0088164  +  .0352655 

6 +  .00959281  +  .06714931 

7 +  .0033255  +  .0199532 

8 +  .00220851  +  .01987641 

9 +  .0015117  +  .0120935 

10 +  .00062171  +  .00683871 

n +  .0007765  +  .0077650 

12 +  .00017721  +  .00230311 

13 +  .0004341  +  .0052094 

14 +  .00006371  +  .00095531 

15 +  .0003296  +  .0046145 

16 -  .00007261  -  .00123471 

17 +  .0001624  +  .0025988 

18 -  .00005501  -  .00104431 

19 +  .0001077  +  .0019383 

20 -  .00005251  -  .00110271 

21 +  .0000694  +  .0013876 

22 -  .00004351  -  .00099981 

23 +  .0000493  +  .0010836 

24 -  .00003831  -  .00095771 

25 +  .0000349  +  .0008364 

26 -  .00003331  -  .00089971 

27 +  .0000251  +  .0006535 

28 -  .00002921  -  .00084561 

29 +  .0000184  +  .0005160 

30 -  .00002591  -  .00080291 
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Figuring  the  stresses  at  z  =  i  vrith  these  values  yields 

C^x  =  *3.073, 

Cy  =  -0.00^3, 

Txy  =  0. 
These  results  agree  qiiite  well  vd.th  those  of  Maunsell  (5)  and 
Isibasi  (3)«  They  also  indicate  that  Frocht  (2)  may  have  been  correct 
in  that  his  results  indicate  that  CJ'x  at  z  =  i  slightly  exceeds 
three. 

There  might  be  some  doubt  as  to  the  use  of  equation  (l.lO)  in  re- 
moving stresses  from  the  send -circle,  since  the  actual  stresses  depend 
on  higher  derivatives  of  y   and  ^  ,     To  satisfy  such  question,  for- 
mulas (3.4a)  and  (3.5a)  will  now  be  applied  directly  to  the  stresses  on 
the  notch. 

Here  it  is  most  convenient  to  change  to  polar  coordinates;  since, 
on  the  notch,  the  radial  stress  is  ^r   and  the  tangential  stress  is  Tf.Q, 
The  relation  used  at  this  point  is  (7,  p.  13S) 

(5.6)  d,-fr..  =  y/'v  yp  -  e"''[Hyr'v  x"] . 

On  the  curve  in  question,  z  =   O  =   e^®,  e^i®  =  O^,     d|.=  0,  and 
(re  °  0»  Then  equation  (5.6)  becomes 

(5.7)  0=  f+y^'-a'[<3-'r+x"]. 

where  l//"  and  X'are  given  by  (5.1).     Substituting  (5.1)   in  (5.7) 
gives 
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-CJ 


0=f  -£na.(^-'^--,|--£na.cj"^' 

OO 

=  |-  +  -f-d^-^[na„+  n(n+i)a^-  (n+2fa„^j  +(n+a)a^,^]<y' 


-h-i 


n-t 


+{a.-a.)  +2(2a^-a^)a" . 

Setting  ^n  ^  '^n  "*"  -^"^n'  ^"^  separating  real  and  imaginary  coeffi- 
cients, we  have 

OO 


00 


-n-i 


r\=i 


The  real  coefficient  part  above  will  be  denoted  by  G{6),   and 
the  imaginary  coefficients  will  be  separated  into  even  and  odd  powers 
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of  ^.  Then 

0=    G(<^)+i^(2n-l)d,„.,a^'^  -iUd) 

n  =  i 

n=l 

-i^[2n(2n+2)d,,-(2n+2Xan+3)d,,^Jcr-'""'  . 

n  =  l 

Now  formulas  (3«4a)  and  (3.5a)  may  be  applied  to  give 

0=  Gto)  *^-ti^-\)'i^.,tj^^  -^d£  -C-'o 

n=i  K«-oo 

*>       CO      3K  .22.    ^  2K 

-^[2n(2n-2)d,, -(2n>2)(2n+3)d,.,e]^-2;^ 


n^  K=-«> 


Satisfying  this  eq\xation  removes  stresses  on  the  notch. 

Next  we  equate  the  coefficients  of  the  various  powers  to  zero,  where 

oo 
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For  the  constant  term,  we  get 


0    =    T   ^E^ndJD^I^     +2  fed, 


n=i 


n=i 


oo 


_     P    .    2    , 


00 


which  is  exactly  eqiiation  (5.2). 

The  coefficient  of  the  second  power  gives 


0  = 


oO 


^^-Ei^ 


ir 


Mr  I 


n-2 


an 


2.  '        111       5       ^        Z_.  (Zn-l)(2n+3) 

=    -P.  +  C    +  2  V  .Jiiinidi^ 

2  '       TrZ_j  (2n-l)(2n+3)      • 
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The  preceeding  equation  is  exactly  equation  (5.4).  Likewise, 
the  other  coefficients  yield  the  same  eqviations.  Therefore,  the 
fonnulas  (3.4a)  and  (3»5a)  may  be  used  either  to  remove  resultant 
force,  or  to  remove  directly  stresses  on  a  boundary  vdth  identical 
results. 

We  next  apply  these  formulas  to  a  problem  which  has  not  been 
solved  heretofore. 


CHAPTER  VI 
THE  SEMI-ELLIPTICAL  NOTCH 

The  Mapping 
One  of  the  advantages  of  the  use  of  con^jleoc  variables  is  that 
regions  which  are  difficult  to  handle  may  be  mapped  conformally  onto 
regions  which  are  easily  handled.  To  solve  the  problem  of  a  semi- 
infinite  sheet  with  a  semi -elliptical  notch  under  uniform  tension, 
the  notch  will  be  mapped  on  the  half  unit  circle.  The  mapping  func- 
tion which  will  be  used  is 

(6.1)  Z  =  C^  4-  -^  , 

where  ^  =£"♦-»  7)   »  c  =  ^-J-^  ,  and  e  =  ^  ^^, 

The  curve  C.~  ^      ^s  the  unit  circle  in  the  ^-plane.  Then 

=  acosQ     ■♦-  bisinG  . 

This  is  an  ellipse  with  major  axis  2a  and  minor  axis  2b  situated  as 
in  Figure  2  on  the  next  page.  The  points  C^  =  ±1  give  z  =±a,  so 
that  the  semi-ellipse  maps  onto  the  semi -circle.  Moreover,  ^  =  +  i 


3h 
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gives  z  =ib,  showing  the  half 
of  the  ellipse  vdth  positive  7| ' 
maps  onto  the  part  of  the  circle       ^   p 
vrith  positive  y.  Then  this  map- 
ping transforms  the  semi-ellipse 
exactly  as  desired. 

The  follovdng  theorem  as  stated  by  Tovmsend  is  useful  in  analyzing 
the  mapping  function  Ci7(l)' 

Given  a  multiple-valued  analytic  function  w  =  f(z)  whose  in- 
verse function  z  =  0(w)  is  single-valued.  If  the  derived 
function  0(w)  has  a  zero  of  order  k,  or  a  pole  of  order 
k  +  2  at  Wq,  then  f(z)  has  a  branch  point  of  order  k  at  the 
corresponding  point  z  .  Also  if  w  =  oo  is  a  pole  of  order 
k  of  0'(w) ,  then  the  corresponding  point  of  the  z-plane  is 
a  branch  point  of  f(z)  of  order  k.  (10,  p.  335) 

Now  z  =  cc  ■•■  ec  is  single -valued,  and  its  derived  function 

(j,7'=C-eL   has  a  pole  of  order  2  at  C  =  0  and  zeros  of  order  one 

at   L=^V®/C  •  Therefore,  by  the  theorem,   ^=  f(z)  has  branch 

points  at  the  z  values  corresponding  to   L  ^iyG/C  .  To  find  these 

points,  the  C  values  are  substituted  in  (^■'.l). 


Z  = 


+  e" 


=  zVce 


-H 


J,  L  ~~\lc    yields  the  value  z 


•Va^  -  b^.  These  are 


Likewise 

the  well-known  values  for  the  foci.  They  may  be  Joined  by  a  branch 
cut,  leaving  the  function  holomoi7)hic  in  the  plane.  These  branch 
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points  lie  outside  the  region  of  the  body  in  question,  and,  therefore, 
cause  no  difficulty  in  this  problem. 

Before  proceeding,  the  x-axis  image  must  be  found  in  the  r-plane. 
To  do  this,  equation  (6.1)  is  solved  for  corresponding  values  in  the 
C  -plane  to  z  =  x  +  iO. 

Equating  real  and  imaginary  coefficients  gives 

e 


c-   fi 


r^' 


71  -    0     . 

The  above  equation  is  satisfied  if  7?  =  0.  Then 

cf +  er'      =    X   , 

and  X  can  take  on  any  value  numerically  greater  than'y /^  so  that  7/  =  0 
obtains  all  values  of  x  outside  the  foci  of  the  ellipse. 

If  C  ~^/c^*+  7?*=0>  ^h®  relations  (6,2)  are  also  satisfied.  Then 

c^%  cTf-  e  =  0   , 

This  evidently  yields  a  circle  of  radius  "A/tt  in  the  ^-plane.  The 
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corresponding  value  of  x  woxild  be 


But 


so  that 


5  '  -  \/  r     •» 


|x|<VcF{2)=V^^rb-^  . 


The  above  values  of  x  fall  between  the  foci  of  the  ellipse  and  out  of 
the  region  of  consideration.  Then  the  part  of  the  x-axis  on  vhich  the 
boundary  lies  is  given  by  the  ^-axis  in  the  ^-plane. 

The  mapping  is  now  ready  to  be  used,  but  it  is  helpful  to  consider, 
before  proceeding,  the  form  that  equation  (l.lO)  takes  in  a  mapped 
region. 

Equation  (l.lO)  in  the  ^ -Plane 
The  formula  (l.lO)  will  be  used  as  before,  although  its  form  will 
be  somewhat  altered  by  the  mapping.  In  the  z-plane  it  takes  the  form 
below  in  which  subscripts  have  been  added  to  simplify  succeeding  work. 


where  H  is  a  complex  constant. 
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Now  z  =ijJ{C)   so  that  (l.lOa)  becomes 


(6.2) 


dz 

f  (^)  +  u(^)v;tidu2  +x^  =  H , 

dz 

w^^ri^)  +  x^  =  H . 


This  form  will  be  useful  after  the  mapping  has  been  applied. 

The  Complex  Potentials  in  the  ^-Plane 
Since  the  tension  at  infinity  and  the  conditions  on  the  x-axis 
are  the  same  as  were  taken  for  the  circular  notch,  the  complex  poten- 
tials take  the  form  in  the  z-plane 

oo 

(2.1b)      y/;  =  ALogE  +-;^Z  +^a^Z""   ^ 

(2.2b)     x.'  =-A  Loc^E  --f-z  +^(na^-aj£-"  . 


r\  =  l 


Since  the  boundary  of  the  problem  is  symmetrical  about  the  y-axis 
and  (l.lOa)  holds,  then  the  logarithmic  term  must  vanish  as  it  did  in 
the  circular  notch.  Then  the  equations  (2.1b)  and  (2.2b)  simplify  to 

p 


After  mapping  the  functions,  we  have 


(6.3) 


Pc   -     ^.   .-n 


^•-"ti  +Eb„^ 


« 

n  J 

nri 
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The  first  degree  terms  arise  from  the  first  degree  terms  of  the 
preceeding  forms*  All  other  positive  powers  vanish  since  in  the  C- 
plane  the  stresses  are  bounded  at  infinity.  The  other  coefficients 
are  difficult  to  ascertain  in  terms  of  the  original  z-plane  coef- 
ficients, so  the  solution  vd.ll  be  carried  out  in  the  r^-plane.  To 
find  the  stresses  in  such  a  solution,  the  C  value  corresponding  to 
a  desired  z  value  is  found,  and  this  is  used  in  the  (^-plane  solution. 

The  Problem 
Now  the  preceeding  work  will  be  applied  to  the  elliptical  notch. 
It  has  been  shown  that  the  ^-axis  corresponds  to  the  x-axis.  Stresses 
will  first  be  removed  from  this  axis.  In  the  application  of  (6.2),  use 
must  be  made  of 


(6.A) 
where 


Mi 


ut] 


c  —  e 

F 


2A^ 


-4- 


Do 

^^ 

1 

^ 

D„ 

= 

2X 

^ 

D„. 

" 

0 

m 

?<--  f 
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Substituting  (6.3)  and  (6J+)  in  (6.2)  and  setting     C""S     si^ea 

h  =  i 

If  we  take  the  eeries  product  by  the  Cauchy  method  and  combine  terms, 
we  have 


oo   ri 


-nEKD„.K]£-". 


where  H  is  a  complex  constant.  Simplifying  yields 

=  na„^b„+D„.(^)-jKap„.jr". 

There  is  no  constant  term  so  that  H  is  zero  as  before.  Equating  the 
coefficient  of  each  power  to  zero  gives 

K  =  l 

If  the  coefficients  conform  to  the  above  relation,  then  the 
^  -axis  will  be  free  of  stress.  Thus  stresses  are  removed  from  the 
^-axis  if 

(6.5)         b„  =  iKaA-K  -a„-D„„(^)  ,  n-.^i,},-. 

K  =  l 


a 


Consideration  must  be  made  of  the  curve  c  =  ^  .  Along  this 
curve,  it  will  be  necessary  to  use 


(6.6) 


vrhere 


u{ 

a 

co'  +  ecJ"' 

(^t) 

_    e<3"'  +  c<J 
C  -  ecJ* 

=  Act-'  +(n-A')cr  +A(IUlo'-^ 

h:0 

< 

'b.  =  A  ,        A=  -§-, 
t3.„=  r{\*X'  ,     n>o, 

..=  0  ,       n>o. 

Setting  ^=^  and  substituting  (6.3)  and  (6.6)  in  (6.2)  leads  to 


0  = 


-f^^+Eb^c," 


r»si 


•-11=0  -"-^,=  1  J       ^_| 


-I 
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Again  taking  the  Cauchy  product  and  simplifying  gives 

0=   ^cJ-^cf-^f:a„cr-"+nb.H^)Bjc^" 


4 


~  h=l  K=l 

and  substituting  equation  (7.4)  into  this  equation  gives 


»sl   K:l 


Into  the  above  expression  is  substituted  a^  =  c^  +  id^^  and 
the  part  with  the  real  coefficients  is  vrritten  first  to  get 


0=  -^^ +[^-^]<^-'  +Ec„ 


d" 


m  K=i 

+£id.<^-"  +  i:i[D-Kdj(D„..-Bj-a>". 

n=i  n=i   K=i 

Let  the  part  with  real  coefficients  be  called  g(c^)  and  let  the  even 
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and  odd  powers  on  the  remaining  part  be  sepeirated.  Then 

oo   2n 


0  =     g«J)  +Zi4''"'V^iE(-kdJ(D„..-B.„.J-d.„]o- 


The  equation  is  now  in  proper  form  for  the  application  of 
formulas  (3.Aa)  and  (3.5a). 


0=    9(«)  +  f  Zd„^^ 


**  °°  2m-i 


(6.7) 


2m -I 

oo     an  ****  2       I 

OO  oo 

^2m 


+Ed„-,x:- 


n=i  w  =  - 

oo      tn-l  oo 


nEKd4a„.,..-B„.,.j+d.„.,]x:,;^;^ , 


+ 


where 


n  =  (       K:| 


From  (6.7)  the  coefficients  may  be  taken  to  form  the  final  equa- 
tions. As  in  the  circxilar  notch,  the  coefficients  of  even  powers  yield 
a  homogeneous  system  in  odd  d's  and  even  c*s.  Moreover,  these  num- 
bered c»s  and  d»s  do  not  occur  in  the  coefficients  of  odd  powers  >jhich 
have  only  even  d»s  and  odd  c*s.  Therefore,  all  even  c»s  and  odd  d*s 
are  set  equal  to  zero,  and  that  much  of  the  system  is  determined. 
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D„-B„  = 


In  simplifying  the  coefficient  eqxiations,  the  following  in- 
formation, taken  from  (6.4)  and  (6.6),  is  needed: 

(  -  X  ,    m  =  o, 

Setting  the  coefficients  of  the  odd  powers  equal  to  zero  gives 
the  following  eqxiations: 


K  =  l 


(6.8) 


-i< 


n:i      K=l 


-K^-^^ 


0=    c^_ 


Ztn  —I 


an 


2    < 

"tT 


E[ 


2(2n)d 


£D- 


(^►,f-(2m-lV 


rr\  -  2,3,4,'--. 
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These  equations  pxxiduce  a  system  of  equations  which  give  the 
desired  result.  To  show  that  they  are  compatible  with  the  semi- 
circle equations,  A  will  be  set  equal  to  zero  and  the  results  in- 
vestigated. Since  A  =  q  +  b  '  A=  0  gives  the  result  when  a  =  b, 
i.e.,  a  circle.  With  this  substitution,  the  equation  for  the  first 
power  coefficient  becomes 


Pc 


0=  ^^^-c^c+l^^^ 

oo      in 


r\-\    KK 


Pc     .     2   , 


Y"    a(2n)dtn    ,  Y^  2ndan  -Un-^ld 


zr-a 


=  ^ 


*f 


n-i  n=A 


This  is  identically  equation  (5.2)  of  the  semi -circle  solution  when 
0=1.  Thus,  this  solution  of  the  elliptical  notch  is  compatible  with 
the  solution  in  Chapter  V. 

The  eqimtions  (6.B)  may  be  written  in  a  form  that  is  better  for 
computation.  The  altered  equations  are 


f* 
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m-i 


(6.3a)  ^.  2 


0=  ^Pcx-'(l-xr-^|aK-i)c_.x'"-^-'(l-xr 

K  =  l 

-+-[2m-a-(am-|)A]Ca„.. 


-'*  0=  -^  +  ^+C. 


2ti 


-''-•*  0=  C 


2f«-l 


These  equations  compose  the  exact  solution  to  the  elliptic  notch.  The 
infinite  series  in  the  coefficient  causes  no  question  since  it  con- 
verges for  all  values  of  A« 

To  find  solutions  of  any  particular  ellipse,  a  finite  system  of 
dimensions  corresponding  to  accuracy  desired  is  taken  from  this  infi- 
nite set.  Suppose  an  oi  x  ot  system  has  been  chosen.  Then  n  is 
Slimmed  from  1  to  -g  in  each  equation  and  k  is  summed  from  1  to 
•y  -  n  in  each  d  coefficient.  The  first  equation  is  taken  from 
the  -1  power  coefficient  of  (6.a),  the  second  from  the  +1  power  coef- 
ficient. The  odd  nimbered  equations  from  3  to  OC-  1  are  taken 
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from  the  -(2m  -  l)  power  coefficients,  and  the  even  nvunbered  equa- 
tions above  two  from  the  (2in  -  l)  power  coefficients. 

When  taken  in  this  way,  the  equations  tend  to  have  the  greatest 
coefficients  on  the  main  diagonal  of  the  coefficient  matrix,  as  did 
the  circular  notch.  This  increases  the  value  of  the  associated  de- 
terminant so  that  the  accuracy  of  the  solxition  is  improved. 


CHAPTER  VII 
CONCLUSION 

Previous  Solution  of  the  Semi-circular  Notch 
The  solution  of  the  semi -circular  notch  problem  by  Maunsell  (5) 
has  been  mentioned.  It  is  interesting  to  consider  the  relationship 
of  his  solution  to  the  one  developed  in  this  paper. 

Maunsell  started  with  the  real  stress  function  (j>   vrhich  placed 
uniform  tension  P  in  the  x-direction.  To  this  he  added  certain 
stress  functions  vfhich  did  not  alter  the  stresses  on  the  x-axis,  but 
did  give  him  something  with  which  to  remove  stresses  on  the  notch. 
The  functions  he  added  were  even  numbered  derivatives  with  respect  to 
X  of  R&  cos  6  and  odd  niombered  derivatives  with  respect  to  x  of 
R0  8in0  .  These  are,  respectively, 

(7.1)     l!!^^^=^,[(2n-()sin(a«*l)0-(£r,>l)sln(2n-l)0]^^>,, 
(^2)      a'""(R0sine)  _  Mi__rcos2n0  -  COS  (Sn +2)91  ,  m>0. 

The  numerical  factor  is  combined  and  called  M  and  M,  for  simplicity. 
Attached  to  each  of  these  was  an  undetermined  coefficient  which  he 
later  evaluated  by  taking  a  system  of  eqioations  from  an  infinite  set. 

Analyzing  this  solution  by  use  of  con^jlex  variables,  we  find 
that  Maunsell' 8  original  stress  function  is 
(7.3)  (f>    =  r{azz  -  Az^}  ,  where  A  =  2. 
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Now  the  functions  that  he  added,  (7.1)  and  (7.2),  can  be  expressed  in 
complex  form  as 

(7.1.)  r{mi[^-^^]},       K>.  , 

(7.2a)       R  {«.[-jW  -  jStt]}  1     K  >  1  . 

Adding  these  fonns  to  (7.3),  attaching  a  real  coefficient  to  each,  and 
allowing  these  coefficients  to  absorb  the  M  factors  gives 

This  separation  suggests  the  use  of  the  bracketed  terms  for  Uf  and 
X  as  complex  potentials.  It  is  now  seen  that  if  we  take  the  first 
bracket  for  Xjf ,   we  have  real  coefficients  on  odd  powers  and  pure  imagi- 
nary coefficients  on  the  even  powers  of  z.  The  more  general  attack 
of  this  paper  shows  that  this  must  be  the  case.  Taking  this  function 
for  Xp   associates  the  bj^  above  with  the  o^^-.   of  Chapter  V.  Also 

(2k  -  l)a  is  associated  with  d  .  These  solutions  show  five  signifi- 
k  2k 

cant  digit  correlation  for  the  first  two  terms  and  decreasing,  but 
good,  correlation  thereafter.  This  is  shown  by  the  san^jles  from  the 
two  solutions  in  the  table  at  the  top  of  the  neixt  page. 
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TABLE  II 
COMPARISON  OF  COEFFICIENTS 

Maunsell  Chapter  V 

h^   =  0.792A661  c^  =  0.79248481 

a-j_  =  .3314167  d2  =  .33141212 

b„  =  .00041  C]^3=  .00043 

lla^  =  .00017  d^2=  •00018 

The  fact  that  the  results  of  this  paper  seem  to  be  more  acciirate 
than  those  of  Maunsell,  judging  from  the  lower  ^y  at  z  =  i  and 
more  coefficients  being  figured,  also  tends  to  lend  evidence  to  the 
photoelastic  results  of  Frocht.  Maunsell  achieved  a  value  of  3*05  P 
for  Cj,  at  z  =  i,  and  predicted  a  value  of  3  P  as  the  exact  an- 
swer (5)»  The  improved  accuracy  of  this  paper  shows  a  value  of 
3.07  which  agrees  with  Frocht 's  prediction  that  the  stress  slightly 
exceeds  3  P  (2).  Isibasi»8  figures  show  3.06  P  (3). 

The  Way  Ahead 
This  work  suggests  that  there  may  be  many  such  expansions  which 
could  be  usefxil  in  solving  multiple  boundary  problems.  At  first  glance, 
these  methods  may  seem  rather  comp]JLcated,  but  they  have  the  advan- 
tage of  giving  equations  in  forms  that  may  be  easily  manipulated.  By 
allowing  the  use  of  conformal  mapping,  those  methods  render  possible 
the  solution  of  problems  not  practically  feasible  by  other  methods. 
Also,  if  other  such  solutions  follow  the  general  trend  of  the  two 
shown  in  this  paper,  they  will  lend  themselves  readily  to  electronic 
calculator  solutions.  With  the  increasing  availability  of  such 
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computers  in  research  and  industry,  solutions  as  these  may  well  soon 
be  used  in  many  practical  applications. 
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